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ABSTRACT 


Polymer insulators tend to fail because of the climatic and environmental 
conditions. The failure occurs when the surface of insulator is contaminated 
by sea salt or cement dust which lead to partial discharge (PD). Leakage 
currents will increase by PD that causes deterioration of insulation. 
To predict the insulation failures, an adaptive neurofuzzy inference system 
(ANFIS) method using initial color detection processes are proposed to 
estimate the leakage currents based on the polymer insulator thermal images 
(infrared signature). In this study, the sodium chloride and kaolin are used as 
pollutants of the polymer insulator according to IEC 60507 standards. 
Then, the insulator is tested in the laboratory using AC high voltage applied 
at 18 kV where the temperature detection is controlled at 26° C and 70% RH 
(relative humidity). The percentage of colors (Red, Yellow, and Blue) from 
the thermal image is measured using the color detection method. Correspond 
to the color percentage, the ANFIS method predicts leakage currents from 


polymer insulators. Furthermore, this system interprets measured data from 
insulators that need to be categorized as Safe, Need Maintenance or Harmful. 
The final application of the system can be a non-contact tool to predict the 
polymer insulators used by technicians in the field. 
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1, INTRODUCTION 

Polymer insulators [1-17] show several advantages over ceramic insulators [18-31] such as low- 
energy surfaces that help to maintain hydrophobic surfaces. However, insulators in electricity transmission 
may not function properly even though the age of isolation 1s still within reasonable limits. This can be 
occured due to insulators work under different climatic and environmental conditions such as near the beach 
or a cement plant. The surface of this polluted insulator [1, 6, 10-12, 14-16, 18-20, 22-31] can create partial 
discharge (PD) [13, 22, 28, 32] and lead to increased leakage current [2-7, 9-16, 21-22, 26-29, 31-32]. 

As leakage current 1s getting increased, it causes a decrease in the quality of the insulator and 
shortens the lifetime of the insulator in the insulation system. If an intensity of the contamination layer on the 
insulating surface increases for a long time, the leakage current may cause heat and even leads to surface 
flashover [1-2, 6-8, 10-12, 14, 20, 25-27, 29]. 
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Different methods have been applied to study leakage current and PD behavior of contaminated 
surfaces [30] on insulators [10, 20-25, 31]. However, a correlation study between leakage current and thermal 
behavior of insulator surfaces has not been found in the literature [6]. In this proposed method, a thermal 
camera [2, 8, 18-19, 28, 32] 1s used to capture an image of thermal surface of the insulator and then the 
thermal image is interpreted using a color detection method. By calculating the percentage of color which is 
used as input in the ANFIS method to predict [1, 4, 26-27, 29] the leakage currents of polymer insulators in a 
non-contact manner. A MatLab program has been developed in the system to classify the condition of the 
insulator into category: Safe, Maintenance and Dangerous state. Based on the category result, a field 
technician can make a decision whether necessary or not to replace the insulator. 


2. RESEARCH METHOD 
2.1. Test Objects 

A 20 kV polymer insulator was used in this work as illustrated in Figure |. The insulator tested was 
hung vertically [10, 12, 14, 19, 21, 23-25, 31] facing the thermal camera in a space (1.2 x 0.5 x 0.7 [m3]) as 
shown in Figure 2. The test voltage was 18 kVrms at a frequency of 50 Hz. Room humidity was at 70% RH 
and temperature was 26 °C. The test was carried out in accordance to IEC 60507; clean fog test procedure. 
Before testing, the insulator surface was cleaned with isopropyl alcohol and rinsed with distilled water to 
remove dirt and grease [27]. Typical beach salt pollution was produced by creating a surface layer of 
contaminants on the insulator. Contaminants were formed from NaCl with 40 grams of kaolin mixed with 1 
liter of water. The NaCl salt concentration gives a density of equivalent salt deposits (ESDD) in mg/cm2 with 
a heavily contaminated value [27]. Liquid potassium sulfate-producing containers were used to maintain 
relative humidity in the fog room. The fog producing container was connected by two tubes to the polymer 
insulator test container. Relative humidity and temperature in the test container were maintained using an 
Arduino microcontroller with a temperature sensor, humidity sensor, two sets of 12 VDC fan for air 
circulation, a 12 VDC motor for potassium sulfate stirrer and a heater. 











(b) 


Figure 1. The polymer insulator of 20 kV under level of (a) clean and (b) contaminated 


2.2. Measurement System 

As shown in Figure 2 that the voltage regulator was connected to the input part of the High Voltage 
Transformer (HVT) which was used as a test voltage regulator applied to the polymer insulator. 
High Voltage Amplifier in the form of (HVT) was used as a high voltage generator which was used as a test 
voltage [26] with a ratio of 1: 2000 Volt (peak to peak). An oscilloscope (Tektronix, DPO 5104) was used to 
display the input voltage via a voltage probe (Tektronix, P6015A) as well as the leakage current through the 
insulator and also the PD activity was recorded by the PD detector. The high voltage of the transformer 
through 25MQ series resistance enters the 4000 pF insulator and capacitor filling. If there was a PD on the 
insulator then the capacitor has removed the charge so that the PD detector would detect the PD occurrences 
and then the PD activity was displayed in the oscilloscope. Voltage probe with a ratio of 1: 1000 V was used 
as a working voltage/test detector applied during testing. Leakage current was calculated based on a voltage 
read by an oscilloscope on a 47 kQ shunt resistor [6, 22, 27] which was in series with the insulator. 
The purpose of using a 47 kQ sampling resistor was that the leakage current signal can be detected properly 
for a high voltage range of 10 kV to 24 kV [21]. Thermal camera from Thermograph FLIR A-600 series was 
used to capture thermal images and then stored and processed to predict leakage currents with ANFIS 
intelligent systems via laptop [27]. 
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Figure 2. A schematic diagram of the experimental setup 





2.3. ANFIS System 

A typical ANFIS architecture where each circle shows a fixed node and each square indicates an 
adaptive node as shown in Figure 3. So the basic system of rules has two rules if-then the Takagi-Sugeno 
type as follows: 

Rule 1: if x is Ai and y is Bi, then fi=pix+qiy+ri, 1=1,2. 


Laver 1 





Figure 1. Architecture of typical ANFIS 


In Figure 3. there are five process layers, namely: 
Layer 1: fuzzification process 
Each node 1 1s an adaptive node and the output is given by 


Ol,i = wAi(x), for i= 1, 2, or (1) 
Ol,i=uBi-2(y), for i = 3,4, 


The membership function for A can be in the form of the appropriate parameterized membership 
function, such as the generalized bell function. 


A(x) = se 


1+|(x-c,)/a, 





(2) 


Layer 2: product layer 
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Each node is a fixed node labeled a that determines the activation strength (firing strength) of each 
rule as a product of all incoming inputs, 


O2,1 = 1 = wAi(x)uBI (y), 1 = 1,2 (3) 
So that 
ol =uwAIl(x) AND uBI (y) (4) 


w2 = wA2(x) AND uB2 (y) 


The output in this layer acts as a function of weight 
Layer 3: normalization layer 

Each N label circle node is non-adaptive which calculates the ratio between firing strength on the 
first rule to the total firing strength of all rules. 


W. 
C= — 
ee es (5) 


Layer 4: defuzzification layer 
The output of each adaptive node in this layer is 


O4,i = wifi = wi(pix+qiy+ri), (6) 


where wi is the output on layer 3 and {pi, qi, ri} are called as consequence parameter. 
Layer 5; total output layer 

The Final Layer, the single node in this layer 1s a fixed node labeled, which calculates the overall 
output as the sum of all incoming signals. 


-0,,.-Y wf - sah 


L 


Overall output (7) 


3. RESULTS AND ANALYSIS 

With a heavily polluted insulator at 18 kV AC applied high voltage for about 1 hour (529 data files), 
it would produce thermal on the surface of the insulator against time in the test container at 70% RH and 
temperature of 26 °C. Thermal colors of the surface of the insulator were grouped into 3 main colors: 
blue (303.4° K <x < 304° K), yellow (304° K< x < 306.5° K) and red (> 306.5” K). 

Figure 4a shows the test results in the form of thermal photos of polymer insulators which are 
proportional to the percentage of thermal colors (Figure 4b) and thermal distribution (Figure 4c). 
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(b) (c) 


Figure 4. Testing of heavy pollutant isolators (a) Thermal Photo Zoom (b) Color Percentage 
(c) Graph of Thermal vs Number of thermal pixels 
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Based on the thermal photograph of the insulator surface, the red color is wider in proportion to the 
thermal rise. This can occur because the ridge can accelerate the greater red color of the surface of the 
insulator or the higher the thermal peak of the polymer insulator as shown in Figure 5. 





(a) (b) (c) — d) 


Figure 5. Thermal photos and 3-dimensional surface images of insulators for leakage current 
(a) small (4.2 mA) (b) medium (14 mA) (c) high (24 mA) (d) very high (47 mA) 


3.1. Color Detection Method 

By using the color detection method with the MATLAB programming language, the thermal image 
of the insulator area was shown in Figure 5a with 175,488 (384 x 457) pixels. The first process of Figure 5a 
as the original image is stored in the jpg extension and named as isolatorsmall.jpg as shown in Figure 6a. 
Then, the color detection process approaches red color with the following equation: if (am (i, j, 1)> 100 && 
im (1, j, 2) <90 && im (1, j, 3) <90) with results as shown in Figure 6b. Furthermore, Figure 6b is stored in 
the png extension, which was named as isolatorsmall.png. After that, it was done to return the results in 
Figure 6b to a binary image where the red pixels were replaced with white and the others (which are not red 
pixels) are replaced with black as shown in Figure 6c with an area of 381.405 (541 x 705) pixels. Binary 
images were needed to convert RGB images into binary images (black = 0 and white = 1). Finally, the binary 
image was converted to a binary image negation as in Figure 6d, resulting in a red area of 3,374 pixels. 





(a) (b) (c) (d) 


Figure 6. Detect the red color of the insulator thermal image (a) the original image (b) the color detection 
results are close to red (c) the results of the binary image (d) the result of the binary image negation 


In the same way as Figure 6, the yellow area was calculated by the equation if (im (, j, 1)> 100 && 
im Gi, j, 2)> 100 && im Gi, j, 3) <90) which equivalent to 15,557 pixels as shown in Figure 7. Meanwhile, 
the color area of the insulator as shown in Figure 8 covers an area of 74,519 pixels with the percentage of 
thermal insulator color as plotted in Figure 9 by using the formula; the area of blue = the area of the insulator 
- (red area + yellow area). 
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Figure 7. Detect yellow color of insulator thermal image (a) original image 
(b) color detection results are approaching yellow (c) binary image results (d) binary image negation results 





(a) (b) (C) (d) 


Figure 8. Color detection of insulator thermal images (a) original image (b) the result of color detection 
approaching the color of the insulator (c) the results of binary images (d) the results of binary image negation 


The color detection method above is based on the Otsu algorithm [19] so that the thermal color 
percentage of the insulator surface is obtained as shown in Figure 9. 
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Figure 9. Thermal percentage 


In the thermal percentage for the four thermal conditions of the surface of the insulator, if the greater 
the percentage of red then on the contrary, the smaller the percentage of blue as shown in Figure 9. 
Percentage of yellow remains, this means that PD activity in Figure 10 1s still large within 1 hour thereby 
resulting in the hotter surface of the insulator thus increasing the leakage current as shown in Figure 11. 
The heat from the surface of the insulator disappears and the dry band has formed because the flow of 
leakage current has caused water vapor to evaporate from the surface of the insulator [14]. If the dry-band 
discharge is continuous, the surface becomes more porous which in itself reduces hydrophobicity and partial 
expansion of the discharge (PD) during some dry-bands that eventually produce flashover [15]. 
With contamination deposition on the surface of the insulator, the more influence on the flow of leakage 
current would occur due to hydrophobic reduction [22]. 


XY Geaph 3 Pict 0 Lal 
mo 





(b) (0) 


Figure 10. PD patterns at leakage current (a) 4.2 mA (b) 14 mA (c) 24 mA (d) 47 mA 
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Figure 11. Classification of leakage current (a) 4.2 mA (b) 14 mA (c) 24 mA (d) 47 mA 


With the increasing amount of leakage current as shown in Figure 11, the smaller the phase 
difference between the test voltage and the leakage current, the more resistive it would be [28] which is 
shown in Figure 12. As long as the surface of the insulator is clean, the leakage current shows the behavior of 
the RC network (the leakage current leads the test voltage by the phase angle of around 70 degrees). 
Similar behavior is also shown during the condition of dry polluted surfaces. Therefore, when the insulator is 
dry, regardless of the pollutant layer, the insulator behaves like a high impedance RC network, which allows 
a very low flow of leakage current, <1 mA [10]. Figure 12a shows the test of the net isolator with a large 
enough phase difference between the test voltage and the leakage current at 71% RH at 26.4 °C with a test 
voltage of 20 kV which produces a leakage current of +1.5 mA. When the insulator is wet and lightly 
polished, a significant decrease in phase angle [14] (about 30 degrees). This clearly shows that the resistance 
of the insulator surface is greatly reduced during wet conditions and the leakage current increases when it is 
resistive [10]. When the insulators are polished at a given high voltage, the conductivity of the surface where 
the electric field changes from capacitive to resistive [16]. Figure 12b shows the test of heavy polluted 
insulators at 70% RH, and temperature of 26 °C with a test voltage of 18 kV and the test time has taken about 
40 minutes so that the surface of the insulator has dried which results in a leakage current of 4.2 mA. 
During short use and long arc formation, and the surface 1s heavily polluted and wet, phase angle values have 
become as low as 5 to 15 degrees [10]. As with the test in Figure 13b, Figure 12c shows a small phase 
difference with a leakage current of 47 mA at the instantaneous time of 6.2 minutes of the initial test. 
The graph of leakage current during the test has decrease with the dry rate of the insulator surface at a fixed 
relative humidity of 70% and temperature of 26 °C as shown in Figure 13. Along with the decreasing leakage 
current which has been identified to be directly related to the decreasing thermal surface of the insulator, 
thus PD activity has also getting smaller at the end of the test. 
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(a) (b) (c) 


Figure 12. Phase classification of the leakage current and reference voltage of the insulator 
(a) clean (b) pollutants at 4.2 mA leakage current (c) pollutants at 47 mA leakage current 





(a) (b) 


Figure 13. Heavy polluted insulator at 70% RH, temperature of 26 °C and 18 kV test voltage 
(a) leakage current (b) PD 
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When the humidity is high, the insulator becomes moist, thereby causing PD to occur on the other 
hand rising the leakage current. In addition, the peak value of the leakage current from the polymer insulator 
tends to increase with the increase in relative humidity and the increased pollutants [23]. 


3.2. Grouping of Leakage Current 
Leakage currents are classified into 3 stages, which are secure 0 + 50 mA (O + 24 minutes), forecast 
50 + 200 mA (24 + 48 minutes) and danger > 200 mA (after 48 minutes). Flashover occurs at the danger 
stage so the contamination flashover must be prevented as soon as possible during the secure stage [29]. 
To avoid flashover, the condition of the insulator detected through thermal images associated with leakage 
currents at the secure stage can be classified into 3 categories of categories in Figure 14, namely; 
a) The condition is secure if the length (number of lengths) < 0.5 and width > 0.04 the surface of the 
insulator is reddish. 
b) Mazntenance is needed if the length (number of lengths) ranges from > 0.5 to < 0.67 and width > 0.04 
the surface of the insulator is reddish. 
c) Danger if the length (number of lengths) ranges from > 0.67 to < 0.83 and the width of > 0.02 (0.25 
mm) the insulator surface is reddish. 


Figure 14. Leakage currents are classified into 3 categories 


The length and width of the reddish insulator surface of Figure 5 based on the color detection 
method are like Table 1. 


Table 1. Data on the Length and Width of the Insulator Surface are Reddish 


Size 


Thermal icone Wide Category 
small 0,14 0,2 Secure, length < 0,5 
medium 0,5 0,14 Secure, length < 0,5 
high 0,7 0,1 Danger, length > 0,67 
Very high 0,75 0,16 Danger, length > 0,67 


Based on Table 1, it can be concluded that when the thermal is small (4.2 mA current leakage) and 
moderate (14 mA leakage current) includes a safe while the thermal category is quite large (24 mA leakage 
current) and large (47 mA leakage current) is considered dangerous. When thermal insulators are categorized 
as dangerous because the leakage current is close to 50 mA [29]. Guided by the classification of leaky 
currents [29], in experiments with high-voltage surface insulators and tests conducted for more than 48 
minutes so that flashover occurs. However, the flashover does not occur because the leakage current is 
getting smaller because the thermal surface of the insulator is also getting smaller. This is because when 
testing the relative humidity (RH) the insulator surface cannot be wetter with an 18kV test voltage. 
Although one of the factors that will make flashovers occur when the insulator surface 1s very polluted. 


3.3. Forecasts of Leakage Current with ANFIS 


By using combined fuzzy and neural network (ANFIS) method, the leakage current of polymer 
insulator can be estimated with the ANFIS input data namely arusbocor.dat as shown in Table 2. 
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Table 2. Percentage of Color Data and Output of Insulator Leakage Current 


No Thermal images (%) Leakage 
Red Yellow Blue current (mA) 

1 0 0 100 0.15 

2 3 17 80 7.99 

3 6 22 72 11.64 
4 15 19 66 17.03 
5 21 21 58 35.50 
6 26 31 43 51.66 


ANHIS has resulted minimum RMSE training = 0.000124677 which requires 27 neurons to reach a 
value close to zero [24] as shown in Figure 15. 
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Figure 15. ANFIS Program (a) RMSE training (b) FIS resulted 


The results obtained are appropriate and the forecasted leakage currents in heavy polluted insulators 
is 23.7491 mA at percentage of red color (17%), yellow color (20%) and blue color (63%) is shown in 
Figure 16. 

In addition, Sugeno fuzzy program has been made by using ANFIS as typed accordingly to 
‘arusbocor_anfis' then MATLAB would create FIS. Furthermore, the thermal percentage input using evalfis 
instruction was tested and presented in Figure 14 as well as its estimated leakage current values. 





>> fis=readfis('arusbocor anfis') >> fis=readfis('arusbocor_anfis') >> fis=readfis('arusbocor_anfis') 
fis = fis = fis = 
struct with fields: struct with fields: struct with fields: 
name: 'arusbocor anfis' name: ‘arusbocor anfis' mame: ‘arusbocor anfis' 
type: 'sugeno' type: "'sugeno' type: 'sugeno' 
andMethod: 'prod' andMethod: 'prod' andMethod: 'prod' 
orMethod: ' probor ' orMethod: ’ probor * orMethod: ' pre »bor' 
defuzzMethod: 'wtaver' defuzzMethod: ‘'wtaver' defuzzMethod: ‘'wtaver' 
impMethod: 'prod' impMethod: 'prod' impMethod: '‘'prod' 
aggMethod: 'sum' aggMethod: 'sum' aggMethod: 'sum' 
input: [13 struct] input: [1*3 struct] input: [13 struct] 
output: [1x1 struct] output: [11 struct] output: [1xl struct] 
rule: [127 struct] rule: [1*27 struct] rule: [(1*=27 struct] 
>> arusbocor=evalfis([{15 19 66], fis) >> arusbocor=evalfis([17 20 63],fis) >> arusbocor=evalfis([0 3 97],fis) 
arusbocor = arusbocor arusbocor = 
17.0301 23.7491 0.1872 


Figure 16. Evalfis instruction (a) leakage current is obtained according to training 
(b) tested in the form of forecasted leakage current and (c) forecasted leakage current of the clean insulator 


4. CONCLUSION 

Pre-processed uses a color detection method for images of thermal polymer insulators that will 
speed up the process of estimating leakage currents in a non-contact manner. Pre-processing produces a 
percentage of red, yellow and blue as input from the ANFIS method. In line with that, it 1s necessary to 
classify leakage currents to determine the classification of insulators in the Secure, Forecast or Danger stages. 
For the classification of the Secure stage (current leakage <50 mA) it is also necessary to classify leakage 
isolator currents in the category of Safe, Maintenance or Danger. Safe stage classification with Maintenance 
and Danger categories is the time that technicians can use to prevent flahover from insulators. This time it 
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can be used by technicians to clean insulators from pollutants. The leakage current forecasting based on the 
thermal image of polymer insulator using ANFIS method was successfully applied and the ANFIS method 
can speed up the recognizing process of the leakage current value with an acceptable accuracy. 
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